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Abstract. The solar tachocline is shown as hydrodynamically stable against nonaxisymmetric 
disturbances if it is true that no cos 4 # term exists in its rotation law. We also show that the 
toroidal field of 200 Gauss amplitude which produces the tachocline in the magnetic theory 
of Riidiger & Kitchatinov (1997) is stable against nonaxisymmetric MHD disturbances - but 
it becomes unstable for rotation periods slightly slower than 25 days. The instability of such 
weak fields lives from the high thermal diffusivity of stellar radiation zones compared with the 
magnetic diffusivity. The growth times, however, result as very long (of order of 10 5 rotation 
times). With estimations of the chemical mixing we find the maximal possible field amplitude to 
be ~500 Gauss in order to explain the observed lithium abundance of the Sun. Dynamos with 
such low field amplitudes should not be relevant for the solar activity cycle. 

With nonlinear simulations of MHD Taylor-Couette flows it is shown that for the rotation- 
dominated magnetic instability the resulting eddy viscosity is only of the order of the molecular 
viscosity. The Schmidt number as the ratio of viscosity and chemical diffusion grows to values 
of ~ 20. For the majority of the stellar physics applications, the magnetic-dominated Tayler 
instability will be quenched by the stellar rotation. 
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1. Introduction 

We ask for the stability of differential rotation in radiative stellar zones under the 
presence of magnetic fields. If the magnetic field is aligned with the rotation axis then 
the answer is simply 'magnetorotational instability' (MRI). If the field is mainly toroidal 
(the rule rather than the exception) the answer is more complicated. Then the Rayleigh 
criterion for stability against axisymmetric perturbations of Taylor-Couette (TC) flow 
with the rotation profile J? = fi(R) reads 

Hence, almost uniform fields or fields with cx l/R (a current-free field) are stabilizing 
the TC flow and no new instability appears. 

More interesting is the question after the stability against nonaxisymmetric perturba- 
tions. Tayler (1973) found the necessary and sufficient condition 

±(RBl)<0 (1.2) 

for stability of an ideal fluid against nonaxisymmetric perturbations. Now almost ho- 
mogenous fields are unstable while the fields with ocl/R are stable. 

We have probed the interaction of such stable toroidal fields with stable flat rotation 
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laws and found, surprisingly, the Azimuthal Magneto-rotational Instability (AMRI) which 
for small magnetic Prandtl number scales with the magnetic Reynolds number Rm of 
the global rotation similar to the standard MRI (Rudiger et al. 2007). 

In the following, as an astrophysical application of these nonaxisymmetric instabilities 
the magnetic theory of Rudiger & Kitchatinov (1997, 2007) of the solar tachocline is 
presented. In the last Section we return to a TC flow under the presence of strong enough 
toroidal fields presenting first results of the eddy viscosity and the turbulent diffusion of 
chemicals for the Tayler instability (TI). 

2. Solar tachocline 

The tachocline is the thin shell between the solar convection zone and the radiative 
interior of the Sun where the rotation pattern dramatically changes. 

ROTATION TOROIDAL FIELD 



Figure 1. The tachocline formation on the basis of a fossil poloidal field of 10~ 4 Gauss confined 
in the radiative solar interior (Rudiger & Kitchatinov 1997). Left: the stationary rotation profile; 
right: the isolines of the resulting toroidal field with its amplitude of 200 Gauss. 

The nonuniform rotation of the solar convection zone - which is due to the interaction 
of the convection with the global rotation - has no counterpart in the solar core but 
the convection zone rotates in the average with the same angular velocity as the interior 
does. The radial coupling is thus large. This phenomenon cannot be explained by viscous 
coupling (the viscosity below the convection zone is by more than 10 orders of magnitude 
smaller) but it can be explained with a weak fossil poloidal field which is confined in the 
solar radiative interior. For the amplitude of this field only values of order mGauss are 
necessary resulting in a tachocline thickness of about 5% of the solar radius (Fig. [IJ . The 
resulting toroidal field amplitude inside the tachocline of about 200 Gauss mainly depends 
on the magnetic Prandtl number (and the rotation velocity) for which Pm = 5 • 10 -3 
has been used in the model (Rudiger & Kitchatinov 1997, 2007). One can estimate the 
resulting toroidal field in terms of the Alfven velocity Va — / ^//i p simply as 

V A = VP^Uv (2.1) 

with Uq the linear velocity of rotation, Uq = Rfi. For Pm ~ 1 the resulting toroidal field 
amplitude would be of order 100 kGauss which is certainly unstable. For Pm ~ 10~ 4 the 
field strength is reduced to only 1 kGauss so that we have carefully to check its stability. 
The magnetic theory holds for two main conditions: i) the field must completely be 
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Figure 2. Neutral hydrodynamic instability lines for the rotation law without cos 4 #-term in 
3D. Only Al and S2 modes are unstable. The Al mode for wavenumber k — > reproduces the 
Watson (1981) result a = 5Q/S1 ~ 0.286. 

confined in the radiation zone and ii) the magnetic Prandtl number must be small enough. 
There are several possibilities to fulfill the first condition (cf. Garaud 2007; Riidiger & 
Kitchatinov 2007) which, however, shall not be discussed in the present paper. 



2.1. Hydrodynamic stability 

To fulfill the second condition the radiative tachocline must hydrodynamically be stable. 
On the first view this should not be a problem. As the differential rotation in latitude 
forms a shear flow its amplitude, 5(7 = i? C q — ^poic7 decides the stability properties. By 
use of a 2D approximation which ignores the radial coordinate Watson (1981) derived for 
ideal fluids the condition SSI/ fi < 0.286 for stability. This rather large value would lead 
to a stable tachocline. The radial velocity components, however, are small but not zero. 
Also the latitudinal profile of the angular velocity is more complicate than the simple 
cos 2 6*-law used by Watson. We have thus to rediscuss the stability of shear flows with 



(2(0) = O 1 1 - o((l - /) cos 2 9 + f cos 4 9)j (2.2) 

where 8Q/ Qq ~ a; f is the contribution of the cos 4 9 term which describes the shape of 
the rotation law in midlatitudes. At the solar surface we have a ~ 0.286 and / = 0.55. 

The radial density gradient forms a 'negative' buoyancy leading to damped oscillations 
with the frequency 



C p dr J 

S is the entropy. The frequency N in the upper solar core is by more than a factor of 
100 larger than the rotation frequency fl. So the equation system 

du ( 1 \ ' 

— + (UV)u + (uV)U = - -Vp + vAu 

dt \p J 

T\jfi + (UV)s + (uV)Sj = C pX Ar (2.4) 

must be solved for the flow perturbation u and the entropy fluctuation s. It is divw = 0; 
s = —Cpp'/p. T' and p' are the fluctuations of the temperature and the density, resp. 
The mean flow U is given by f|2 . 2[) . The equations are solved with a Fourier expansion 
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exp(i(fcr + m<f> — tot)) in the short-wave approximation kr 3> 1 with m as the azimuthal 
wave number. As usual, in latitude a series expansion after Legendre polynomials is used. 

In both latitude and longitude the modes are global. The parameter including the 
density stratification is 

*=iS^ (2 - 5) 
so that A — * oo reproduces the 2D approximation by Watson (1981) and Cally (2001). 
They showed that only nonaxisymmetric modes with m = 1 can be unstable and the 
same is true in the present 3D approximation. The modes with ug antisymmetric with 
respect to the equator are marked with Am and the modes with ug symmetric with 
respect to the equator are marked with Sm . 

Let us start with the rotation law (/ = 0). The Prandtl number is fixed as 

Pr = - = 2 • 10" 6 , (2.6) 

X 

where v and \ are viscosity and thermal conductivity. 

The main result is given in Fig. [2] which shows the neutral-stability lines for various 
a = 5Q / Q. Only Al and S2 are obtained as unstable (SI is stable!). For k — > the Watson 
result a = 0.286 is reproduced. With radial stratification, however, this critical value is 
reduced to a = 0.21. For ideal fluids Cally (2003) found instability for a — 0.24 which 
also fits our result. Hence, for a < 0.21 the solar tachocline remains hydrodynamically 
stable. There is thus no shear-induced turbulence. Note also that for N — > (mimicking 
the convection zone) no instability exists. 
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Figure 3. Hydrodynamic-stability map for various values / for the power of the cos 4 (9-term 
in the rotation law (|2.2[) . This term destabilizes the SI mode. For high / already rather small 
shear values become unstable. 

The calculations have been repeated with the cos 4 9-term included in the rotation law. 
Then the SI mode becomes dominant and reduces the critical shear. For / = 0.5 the 
maximum shear for stability results to a — 0.16 (Fig. [3]). 

This is a rather small value. Both our modifications of the Watson approach are desta- 
bilizing the differential rotation. If the shape of the surface rotation with / = 0.5 would 
be conserved through the convection zone and the tachocline then only 50% of the surface 
shear can indeed produce hydrodynamic turbulence in the tachocline. If the shear is not 
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conserved (/ = 0) then the minimum shear for turbulence is 0.21, large enough to ensure 
stability. 

A dramatic stabilization, however, of the latitudinal shear results from the inclusion of 
the tachocline rotation law as observed with its large radial gradients into the calcula- 
tions. With a 3D code without stratification one can show that in this case all rotation 
laws with a < 0.5 should be stable (Arlt et al. 2005). 

In order to decide the stability problem of the solar tachocline informations are needed 
about the exact shape of the rotation law there. Charbonneau et al. (1999) analyzed the 
helioseismic data and found a ^ 0.15 and / — 0. After our analysis (and also after theirs) 
such a rotation law is hydrodynamically stable. Obviously, we have to know in detail the 
space dependence (and time dependence) of the internal solar rotation. 

2.2. MED stability 

We consider the tachocline as hydrodynamically stable. The magnetic Prandtl number 
is thus the microscopic one for which we use Pm = 0.005 as a characteristic value. For 
the instability of toroidal fields the ratio of Pm and Pr, the Roberts number 

q=~ (2.7) 
t] 

plays a basic role. For the Sun the typical value 2500 is used. The interaction of rotation 
and toroidal magnetic fields (of a simplified structure) will be demonstrated with these 
parameters. The general result is that for small q the field must be strong to become 
unstable while for q > 1 much weaker fields become unstable. However, the growth 
times (in units of the rotation period) for the weak-fields are much longer (by orders of 
magnitudes) than the growth times of strong fields (of order of the Alfven period) . 

The left plot of Figs. 0] concerns a field with only one belt which peaks at the equator. 
The Alfven frequency J?a is defined by 

B<f, = rsinOyJfiopttA, (2.8) 

and is considered as constant (see Cally 2003). The global rotation is assumed as rigid. 
For this case Fig. 2] (left) reveals the fields with i?A > J? as always unstable with growth 
rates >J?a- Toroidal fields with amplitudes Va > U oq , i.e. ~10 5 Gauss (for the Sun) 
cannot stably exist in the radiative solar core. Even weaker fields with /?a ~ 0.005/2 can 
be unstable but only for very high heat-conductivity. For q = the weak fields are stable. 
For q> 1 they are indeed unstable but with growth rates smaller than 10 -4 . The growth 
time Tgrowth in units of the rotation period is 

Tgrowth = J_ , 2 g) 

T rot 27T7 

with the normalized growth rate 7 = 3(oj)/J?o- The maximum growth time is then > 1000 
rotation periods (but much shorter than the diffusion time). Even for very large q there 
is a magnetic limit below the magnetic field is always stable. From Fig. [4] (one or two 
belts) we find the minimum field is i?A — 0.005/?. For the Sun, therefore, the maximum 
stable field in the model is ^600 Gauss. 

The model is insofar correct as the microscopic diffusivity values (viscosity, magnetic 
diffusivity, heat-conductivity) have their real amplitudes (for ideal fluids, see Cally 2003). 
The model, however, is insofar not correct as the radial profiles of the fields are assumed as 
nearly uniform. Arlt et al. (2007) work with a 3D code without buoyancy (q — > 00, Pm = 
0.01) for toroidal field belts with strong radial gradients and find instability for weak 
fields with amplitudes of order 10 Gauss. 
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Figure 4. Normalized growth rates for the magnetic instability in rigidly rotating stars. The 
curves are marked with their values of the Roberts number (|2.7[) . Left: one magnetic belt peaking 
at the equator; right: two magnetic belts with zero-field at the equator. There are two instability 
domains: for strong magnetic fields ((2a > fi) the growth rates are high (~ Ha/O) and for weak 
magnetic fields (J?a < O) they are very small (~ (Qa/ Q) 2 )- The latter domain only exists for 
q> 1 (perfect heat conduction). 
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Figure 5. A star spinning down moves from the top (stable area) to the bottom (unstable area). 
It is shown that a toroidal magnetic field of 200 Gauss is stable for the solar rotation period but 
the stability is lost for slower rotation. 

Not surprisingly, for two belts with equatorial antisymmetry (the field vanishes at the 
equator) there are some differences to the one-belt model, but the maximal stable field 
amplitudes are always of the same order (Fig. 01 right). 

2.3. Effect of stellar spin-down 

The question arises whether a solar-type star is always able to form a tachocline. The older 
the star the slower its rotation. Hence the rotational quenching of the Tayler instability 
becomes weaker and weaker for older stars so that the instability becomes more efficient. 
By its spin-down the star moves to the right along the abscissa of both the Figs. |U The 
(slow) magnetic decay goes in the opposite direction; this effect is still neglected. We 
assume that the total amount of the latitudinal differential rotation remains constant 
during the star's spin-down. This is a well-established assumption (see Kitchatinov & 
Rudiger 1999; Kiiker & Stix 2001). 
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Figure [5] shows the results. We have computed the normalized growth rates of the 
magnetic instability of rotating stars with a toroidal field of 200 Gauss and a differential 
rotation of 5fi ~ 0.06 day -1 (the solar value). The rotation period is normalized with 25 
days in Fig. [5] so that the horizontal yellow line represents the Sun. In the upper part of 
the plot the 200 Gauss are stable while in the lower part of the plot they are unstable. 
Obviously, the Sun lies in the stable area but very close to the instability limit. We are 
thus tempted to predict that G2 stars older than the Sun (or better: of slower rotation) 
should not have a tachocline. When the toroidal field becomes unstable then the resulting 
turbulence is able to destroy the tachocline rather fast. 

2.4. Chemical mixing 

The flow pattern of the magnetic instability also mixes passive scalars like temperature 
and chemical concentrations. The instability, therefore, could be relevant for the so- 
called lithium problem. In order to explain the observed lithium concentration at the 
solar surface one needs a turbulent mixing beneath the convection zone which enhances 
the microscopic value of the diffusion coefficient of 30 cm 2 /s by (say) two orders of 
magnitude. Note the smallness of this quantity; only a very mild turbulence can provide 
such a small value of the diffusion coefficient 

D T ~ (u' 2 )r corr . (2.10) 

This relation is used here as a rough estimate, a quasilinear theory of turbulent mixing 
has been established by Riidiger & Pipin (2001) also for rotating turbulences. For a 
correlation time of the order of the rotation period the desired mixing velocity is only 1 
cm/s. 

One can estimate the characteristic time by r corr ~ l 2 /Dt with I as the radial scale 
of the instability and Dt ~ 10 4 cm 2 /s. For the radial scale the value 1000 km has been 
found by Kitchatinov & Riidiger (2008). With this value it follows r corr ~ 10 12 s which 
corresponds to a very small normalized growth rate of T rot /r ~ 10~ 6 . The resulting 
toroidal field which fulfills this condition is smaller than 600 Gauss (Fig. [5]). Stronger 
fields would produce a too strong mixing which would lead to much smaller values for 
the lithium abundance in the solar convection zone than observed. 

Our result in connection with the observed lithium values also excludes the possi- 
bility that some dynamo works in the upper part of the solar radiative core. If such 
a ('Tayler-Spruit') dynamo existed then the resulting toroidal fields with less than 600 
Gauss are much too weak to influence the magnetic activity of the Sun with magnetic 
fields exceeding 10 kGauss. 

3. Tayler instability in Taylor-Couette systems 

To simplify matters we consider the pinch-type instability in a Taylor-Couette system 
filled with a conducting fluid which is nonstratified in axial direction. The stationary 
rotation law between the cylinders is fl = a + bjR 2 where a and b are given by the fixed 
rotation rates of the cylinders. In a similar way the stationary toroidal field results as 
Bfj, = AR + B/R. In the following we have fixed the values at the cylinders to il out — 
0.5/?i n and mostly [Ib = B out /B m = 1 is used. The outer cylinder radius is fixed to 2R m . 
Reynolds number Rm and Hartmann number Ha are defined as 

Rm= ^il ; Ha = BinRin , (3.1) 

the magnetic Prandtl number here is always put to unity. 
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Figure 6. Growth rates in units of Q for the two-belts model and for a fixed radial scale of 
1000 km. The right-hand scale gives the estimated values for the diffusion coefficient for passive 
scalars and the uppermost scale gives the magnetic field amplitude in kGauss. 



A detailed description of the used nonlinear MHD code for incompressible fluids is 
given by Gellert et al. (2007). In the vertical direction periodic boundary conditions are 
used to avoid endplate problems. In this approximation the endplates rotate with the 
same rotation law as the fluid does. The height of the virtual container is assumed as 
6D with D the gap width between the cylinders. The cylinders are considered as perfect 
conductors. The code was first tested for the nonaxisymmetric AMRI which appears 
if stable rotation laws and stable toroidal fields (current-free, fis = 0.5) are combined 
(Rudiger et al. 2007). Figure [7] (left) shows the instability domain (solid) which for given 
(supercritical) magnetic field always lies between a lower Reynolds number and an upper 
Reynolds number. For too slow rotation the nonaxisymmetric modes are not yet excited, 
but for too fast rotation the nonaxisymmetric instability modes are destroyed so that 
the field becomes stable again. 



AMRI {instability of m=1 ) for Pm=1 , (i n =0.5 -| 
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Figure 7. Pm = = 0.5 (current-free). Left: stability diagram from nonlinear simulations. 

Dots: stability, stars: instability; the solid line results from the linear theory. Right: m-spectrum 
for Re = 250, Ha = 110. The m — 1 mode contains 69% of the total magnetic energy. 

Between the limiting Reynolds numbers the instability is no longer monochrome but 
also other modes than m = 1 are nonlinearly excited. The m — 1 mode often contains 
the majority of the total magnetic energy (Fig. right). There are also cases, however, 
where m = and m = 1 contain nearly the same amount of magnetic energy. Note that 
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the nonlinear effects can provide remarkable portions of the energy of the instability in 
form of axisymmctric rolls although the basic instability is a nonaxisymmetric one. 




Figure 8. Left: The instability map for /j,n = 0.5, (ib = 1, Pm = 1. The hatched area is 
stable. Right: The growth rate (normalized with the rotation rate of the inner cylinder) along 
the horizontal line at Re = 500. The instability in the upper domain (AMRI) grows slowly while 
in the lower domain (TI) it grows much faster. 



Now with fis = 1 positive axial currents between the cylinders are allowed. Again there 
are two different instability domains for the nonaxisymmetric modes with m = 1. They 
are separated by a stable domain (Fig. [51 left). The upper one is for fast rotation and weak 
holds (Re > Ha) and the lower one is for strong fields and slow rotation (Ha > Re). The 
growth rates in both domains are very different: The weak-held instability is slow and 
the strong-held instability is very fast (Fig. right). The rotation-dominated instability 
disappears for rigid rotation while the magnetic-dominated instability even exists without 
rotation (Re = 0). The latter one is the Tayler instability (TI) under the stabilizing 
influence of the basic stellar rotation (Pitts & Tayler 1985). The rotation-dominated 
('upper') instability appears to be the AMRI which also exists under the modifying 
influence of weak axial currents in the fluid. Note again that i) too fast rotation Anally 
stops both the instabilities, and ii) its growth rate is very small. One must also stress 
that the presented results only concern the most simple case of Pm = 1. For very strong 
fields the stable domain between AMRI and TI disappears. 

The condition for the existence of TI as given in Fig. [8J (left) is > ^/JI^pRQ while 
the condition for AMRI is RQ > B^/^nop. Both the instabilities, however, only exist if 
the rotation is not too fast. Nonaxisymmetric modes are always stabilized by sufficiently 
fast rotation. 

Ap stars (with 10 kGauss and a rotation period of several days) and neutron stars (with 
10 12 Gauss and a rotation period of 10 ms) are rotation-dominated. Their instabilities 
are not of the Tayler-type. In the following we have thus considered the AMRI in more 
detail. For given Ha (= 500) the eddy viscosity, the diffusion coefficient and the 
Schmidt number 

Sc = w T (3 - 2) 

are computed. The eddy viscosity vt is the ratio of the angular momentum transport 
by Reynolds stress and Maxwell stress and the differential rotation. We find a vt of the 
order of the microscopic value (Fig. [HI left). The maximum exists as the instability - 
as mentioned - disappears for too fast rotation. The averaging procedure concerns the 
whole container so that the values in Fig. [9] are lower limits. 
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Figure 9. Left: the eddy viscosity in units of the microscopic viscosity for /iq = 0.5, /is = 1, 
Pm = 1 and Ha = 500. It grows for faster rotation. Right: the Schmidt number (|3.2|l . 

So far the diffusion coefficient for chemicals could only be estimated by £>t — (u^) / fi. 
It proves to be much smaller than the viscosity. Brott et al. (2008) have shown that for 
too strong mixing the stellar evolution is massively affected. A better theory must solve 
the diffusion equation. 

Accepting this approximation the resulting Schmidt number (|3.2p reaches values of 
20 ... 30 (Fig. [9l right). Obviously, the angular momentum is mainly transported by the 
Maxwell stress while the diffusion of passive scalars is due to only the Reynolds stress 
which is much smaller. Quite similar results have been obtained by Carballido et al. 
(2005) and Johansen et al. (2006) for the Schmidt number of the standard MRI. Maeder 
& Meynet (2005) for a hot star with 15 solar masses and for 20 kGauss find much 
higher values of the Schmidt number (10 6 ). Also Heger et al. (2005) work with magnetic 
amplitudes of 10 kGauss for which J?a < & hence the growth rates are small. 
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